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ABSTRACT 

Wc present a novel approach of identifying the Milky Way (MW) and Andromeda 
(M31) progenitors that could be visible as LAEs at z ~ 6: we couple a snapshot 
from the Constrained Local UniversE Simulations (CLUES) project, that successfully 
reproduces the MW and M31 galaxies situated in their correct environment, to a 
Lyman Alpha Emitter (LAE) model. Exploring intergalactic medium (IGM) ionization 
states ranging from an almost neutral to a fully ionized one, we find that including 
(excluding) the effects of clustered sources the first local group progenitor appears as 
a LAE for a neutral hydrogen fraction xm = 0.4 {xhi = 0.1). This number increases 
to 5 progenitors each of the MW and M31 being visible as LAEs for xhi = 10^^; 
the contribution from clustered sources is crucial in making many of the progenitors 
visible in the Lya, for all the ionization states considered. The stellar mass of the local 
group LAEs ranges between 1O^'^~*M0, the dust mass is between 1O*'^~^'^M0 and 
the color excess E{B — V) = 0.03 — 0.048. We find that the number density of these 
LAEs are higher than that of general field LAEs (observed in cosmological volumes) 
by about two (one) orders of magnitude for xhi = 10~^ (0-4). Detections of such high 
LAE number densities at z '^ 6 would be a clear signature of an over-dense region 
that could evolve and resemble the local group volume at z ~ 0. 

Key words: cosmology: theory, galaxies: individual, galaxies: high redshift; galaxies: 
intergalactic medium 
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1 INTRODUCTION 

The search for high-redshift LAEs has advanced tremen- 
dously since the days of their first successful detections, e.g. 
Lowenthal et al. (1991). Advances in instrument sensitivity, 
refined selection techniques and specific LAE spectral sig- 
natures have led to the confirmed detections of hundreds 
of LAEs in a wide redshift range, between z ~ 2.2 — 6.6 
(e.g. Cowie & Hu 1998; Steidel et al. 2000; Malhotra et al. 
2005; Shimasaku et al. 2006; Kashikawa et al. 2006; Hu et 
al. 2010). 

Over the past few years, LAEs have been used exten- 
sively as probes of reionization, high redshift galaxy evolu- 
tion and the dust content of early galaxies (Dijkstra et al. 
2007; Kobayashi et al. 2007; Dayal et al. 2009; Finkelstein 
et al. 2009; Dayal, Ferrara & Saro 2010; Dayal, Maselli & 
Ferrara 2011). However, scant effort has been devoted to 
establishing a link between such high-redshift observations 
{z ~ 6) to those a.t z = 0, especially in the context of the 
local group. 

The first such attempt was made by Salvadori, Dayal & 
Ferrara (2010), who coupled the semi-analytic code, GAMETE 
(Salvadori, Schneider & Ferrara 2007; Salvadori & Ferrara 
2009), which successfully reproduces most of the observed 
MW and dwarf satellite properties at z = 0, to a LAE model 



developed in Dayal et al. (2009) and Dayal, Ferrara & Saro 
(2010). Though a powerful statistical tool, the semi-analytic 
nature of GAMETE meant that the effects of clustered sources 
on the visibility of LAEs could not be properly considered; to 
alleviate such a problem to some extent, all the calculations 
presented therein were carried out assuming the IGM to be 
completely ionized. 

In this work, we extend the calculations presented in 
Salvadori, Dayal & Ferrara (2010). We use a z ~ 6 snapshot 
from the CLUES runs, which are state of the art constrained 
simulations that build up one of the possible merger histo- 
ries of the local group, including the MW and M31 galax- 
ies. This snapshot is coupled to a LAE model presented in 
Dayal, Ferrara & Saro (2010) and Dayal, Maselli & Ferrara 
(2011), that reproduces a number of data sets accumulated 
for high-z LAEs. The main advantage of using these simula- 
tions lies in the fact that they provide information on both 
the physical properties and the spatial positions of all the 
progenitors of the local group at any given redshift. This 
allows us to calculate the effects of clustered sources on the 
visibility of such progenitors in the Lya, for IGM ionization 
states ranging from a completely neutral (xhi = 0.99) to a 
fully ionized (xhi ~ 10^^) one. 

The main questions we aim to answer with such an ap- 
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proach are: (a) how reionized does the volume containing 
the local group have to be at z ~ 6 for any of its progenitors 
to be visible as a LAE?, (b) where do such LAEs lie with 
respect to the underlying dark matter (DM) distribution?, 
and (c) what are the physical properties of the progenitors 
that are visible as LAEs, and how do they compare to those 
of field LAEs observed in cosmological searches at z ~ 6, 
e.g. by Shimasaku et al. (2006). 



2 THE CLUES SIMULATIONS 

For the calculations presented in this Letter, we use 
the CLUES simulations (http://www.clues-project.org). We 
briefly summarize the simulations here and the interested 
reader is referred to Libeskind et al. (2010), Knebe et al. 
(2010), Libeskind et al. (2011) and Knebe et al. (2011) for 
complete details on the central galaxies, their discs and the 
properties of their z — substructure population, as well 
as details regarding the simulation initial conditions, gas- 
dynamics and star formation. 

The simulations used in this study were run with the 
PMTree-SPH MPI code GADGET2 (Springel 2005) in a cos- 
mological box of size 64/i~^ comoving Mpc (cMpc). The 
runs used standard ACDM initial conditions and WMAP3 
parameters (Spergel et al. 2007) such that Qrn ~ 0.24, 
fit = 0.042, Qa = 0.76, h = 0.73, as = 0.73 and n = 0.95. 
The initial conditions were constrained using observations 
of peculiar velocities and the positions of objects in the lo- 
cal volume (Hoffman & Ribak 1991). Since this method only 
constrains large (i.e. linear) scale structures, a local group 
(defined in terms of the mass, relative distance and number 
of members) was selected from a number of constrained low 
resolution simulations. 

The local group, of volume 2/i~^cMpc at z = 0, was 
re-simulated at a much higher resolution using the prescrip- 
tions given by Klypin et al. (2001), to produce two objects 
resembling the MW and M31, each of which contains about 
10*^ particles within their virial radii. Outside of this re- 
gion, the simulation box was populated with lower resolution 
(i.e higher mass) particles to mimic the correct surround- 
ing environment. In the re-simulated region, the DM and 
gas particle masses were taken to be 2.1 x lO^/i~^M0 and 
4.4 X lO''/i~^M0, respectively. The initial mass function used 
is Salpeter between 0.1 — IOOMq. We used the star formation 
rules of GADGET2 with minor modifications (see Libeskind et 
al. 2010 for details); each gas particle was allowed to undergo 
two episodes of star formation, each time spawning a star 
particle of half the original mass (i.e. 2.2 x 10** H^^Mq). 
After the simulation was run to z = 0, we identified cos- 
mological structures using the MPI-|-OpenMP Amiga halo 
finder (AHF; KnoUmann & Knebe 2009) in each snapshot. 

In order to find the local group LAEs at z ~ 6, we 
identify all the particles within the virial radius of MW and 
M31 at z = 0. These particles are then followed back in time 
and located at z ~ 6, in the simulation box. If bound to a 
structure, we identify the corresponding halo: each of these 
halos was then considered a progenitor of the local group 
at z '-^ 6 and used in the LAE calculations, as explained in 
what follows. 



3 LYMAN ALPHA VISIBILITY 

We start the calculations by obtaining the mass, metallic- 
ity and redshift of formation of each star particle in each of 
the MW/M31 progenitors, in addition to the global physical 
properties (total halo/stellar/gas mass) of each progenitor; 
the age for each star particle is calculated as the time dif- 
ference between its formation redshift, Zform and z ~ 6. We 
find 234 and 240 star forming progenitor haloes of the MW 
and M31 respectively. Considering that stars form in a burst 
after which they evolve passively, we use the population syn- 
thesis code STARBURST99 (Leitherer et al. 1999) to obtain 
the intrinsic spectrum for each star particle depending on 
its age, mass and metallicity; the latter two values refer to 
the mass and metallicity of the star particle at the time of 
its formation. The total intrinsic Lya/continuum luminos- 
ity (1216 and 1375 A in the galaxy rest frame respectively) 
produced by stellar sources for each progenitor is then the 
sum of the Lya/continuum luminosity produced by all its 
star particles. The intrinsic Lya luminosity can be trans- 
lated into the observed luminosity such that La = L]^^ faTa, 
while the observed continuum luminosity, Lc is expressed as 
Lc — U^^ fc- Here, fa {fc) are the fractions of Lya (contin- 
uum) photons escaping the galactic environment and Tq is 
the fraction of the Lya photons that are transmitted through 
the IGM. 

We start by summarizing the calculation of fc'. for each 
progenitor, the dust enrichment is calculated assuming Type 
II supernovae (SNII) to be the primary dust factories (To- 
dini & Ferrara 2001). We assume that each SNII produces 
about O.5M0 of dust (Todini & Ferrara 2001), each SNII 
destroys dust with an efficiency of about 40% the region it 
shocks to speeds ^ 100 kms~^ (e.g. Scab & Shull 1983), a 
homogeneous mixture of gas and dust is astrated into star 
formation, and that dust is lost in SNII powered outfiows. 
We assume a slab-like dust distribution such that the dust 
distribution scale, r^ ~ 0.5rg, where the gas distribution ra- 
dius is calculated as rg = 4.5Ar2oo; the spin parameter is 
taken to be A = 0.05 (Ferrara, Pettini & Shchekinov 2000) 
and r2oo is the virial radius. We then use fa = 1-5/c, as 
inferred for LAEs at z « 6 (see Dayal, Ferrara & Saro 2010 
for complete details of this calculation) . 

Lya photons are further attenuated by the 
H I present in the IGM and only a fraction < Tq < 1 
reach the observer. Since the IGM ionization state 
is largely unconstrained at z ~ 6, we explore 
14 different values ranging from an almost neu- 
tral to completely ionized IGM such that xhi = 
0.99, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 10-^ 10-^ 10"'', 
10"''. For each such initial value of xhi (i.e. the value before 
the effects of the ionization regions built by each progenitor 
are considered), we start by computing the radius of the 
H 11 region each progenitor ionizes around itself; we use a 
H I ionizing photon escape fraction of 2%, following the 
results obtained by Gnedin et al. (2008). However, in reality, 
multiple galaxies generally contribute ionizing photons to 
the same ionized region due to source clustering, which is 
then characterized by an effective radius. Then, within the 
effective ionized region each galaxy is embedded in, the 
total photoionization rate seen by it includes the direct 
radiation from the galaxy itself, from the galaxies clustered 
around it, and from the ultraviolet background. This 
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Figure 1. Sky projections of the progenitors of the MW (upper panel) and M31 (lower panel) at 2 ~ 6 in a radius of 2/i~^cMpc. 
Points and circles show the Lya luminosity of the progenitors of the local group, and the LAE^g respectively, for the initial XHI values 
marked over each projection. The luminosity is color-coded according to the bar shown between 10*^~*^ergs~^ and all progenitors with 
Lc < 10*^ergs~^ are shown with yellow points. While the LAE^g are detectable with current instrument capabilities, all progenitors 
with La > lO'^^orgs"^ are within the sensitivity of future instruments, e.g. the Multi Unit Spectroscopic Explorer (MUSE). The DM 
density field in the 2/i~^cMpc box is shown by contours such that the density is highest in the darkest shaded regions, and reaches the 
average cosmological value for the regions colored white. 



procedure is carried out for each galaxy in the simulated 
volume. Assuming photoionization equilibrium within the 
effective ionized region of each galaxy and forcing xhi to 
attain the assigned global value at the edge of this region, 
we use the Voigt profile to calculate the optical depth, and 
hence Ta for Lya photons along the line of sight. Complete 
details of this calculation can be found in Dayal & Ferrara 
(2011). 

Progenitors of the local group are then identified as 
LAEs based on the currently used observational criterion: 
La ^ 10*^ergs~^ and the observed equivalent width, 
La/Lc ^ 20 A. Progenitors visible as LAEs for any IGM 
ionization state are then referred to as LAE_lg- 



4 RESULTS 

We now answer the first question we posed: how reionized 
does the local group volume have to be at 2 ~ 6 for any of its 
progenitors to be visible as LAEs? We find that including the 
effects of clustered sources, the first LAElg (which is a pro- 
genitor of the MW) appears when the average H i fraction 
is Xhi ~ 0.4; if clustering is not considered, the first LAElg 
appears as late as xhi = 0.1. Progressively lower values of 
Xhi lead to more progenitors becoming visible in the Lya 
as expected: 1 progenitor each of the MW/M31 is visible 
for Xhi = 0.3, 1 (4) progenitors of the MW (M31) are vis- 
ible as LAElg for xhi = 0.1, while the number rises to 5 
each for the MW and M31 for xhi ~ 10~^, as seen from 
the left-most to right-most projections for the MW (M31) 
in the top (bottom) panels of Fig. 1. These results can be 



explained as follows: for a given stellar mass (or star for- 
mation rate) and age, the size of the H 11 region that any 
progenitor can ionize around itself increases with decreasing 
values of xhi', a lower value of xhi in the IGM also results 
in a larger number of clusterers. Both these effects lead to 
a lower value of xhi at each point within the H 11 region, 
which itself is now larger. As a result, the Lya photons face 
a lower optical depth at each point within this ionized region 
and are more red-shifted when they reach the edge of such 
a region; these two effects then boost Ta- As expected, all 
the LAEs become progressively brighter in the Lya as the 
value of Xhi decreases, as seen from a comparison of the 
observed Lya luminosity plotted (panels left to right), for 
decreasing xhi for both the MW and M31, as shown in Fig. 
1. As for the distribution of the LAElg with respect to the 
underlying DM density field, we find that for any value of 
Xhi, LAEs lie in dense DM filaments where the halo masses 
are large enough so that the gas in them can cool and form 
stars, as shown in Fig. 1. 

As for the physical properties of the LAElg, we start 
with the total stellar mass, which is possibly the best con- 
strained physical quantity, when FIR data is available (e.g. 
from Spitzer at 3.6, 4.5/im). The LAElg are amongst the 
most massive progenitors and have stellar mass M* ~ 
1O^'^~*M0; for progenitors with lower masses, a combi- 
nation of dust and IGM transmission reduces La below 
10*^ergs~^. As expected in a hierarchical structure forma- 
tion scenario, LAElg with the largest stellar mass have 
taken the longest to assemble and have stellar mass weighted 
ages f , ~ 150 — 250 Myr as shown in panel (a) of Fig. 2, im- 
plying that many of these progenitors started forming stars 
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Figure 2. Summary of the physical properties of the local group progenitors at z ~ 6. In each Panel, black and red points represent 
the progenitors of the MW and M31 obtained from the CLUES runs. The LAE^q of the MW and M31 are shown by colored circular 
and triangular symbols, respectively. The different colors show the progenitors that become visible as LAEs for different initial IGM 
ionization states: Xm = 0-4 (red), 0.3 (blue), 0.2 (green), 10^^ (magneta) and 10~^ (cyan). As a function of the total stellar mass 
(M«), the panels represent (a) the stellar mass weighted age (t*), (b) the mass weighted stellar metallicity (.Z*), (c) the total dust mass 
(l^dust), (d) the color excess E{B — V), (e) fraction of Lya luminosity transmitted through the IGM (Tq) and (f) the observed Lya 
luminosity (La). The last two values are plotted for XHI ~ 10~^; these then represent the upper bounds on the quantities obtained from 
our model. 



at z ^ 7.6. However, we stress that it is the young stars that 
formed over the last ~ 30 Myr that dominate the stellar lu- 
minosity in the entire wavelength range from the rest-frame 
UV to the near Infrared. It is interesting to see that it is 
not necessarily the progenitors with the largest stellar mass 
that become visible for the largest xhi since the visibil- 
ity sensitively depends on the contribution due to clustered 
sources; such contribution does not correlate with the pro- 
genitor mass, as is explained later. 

As star formation is the main source of metals, progeni- 
tors that are able to build up the largest stellar mass are also 
the most metal rich, as seen from panel (b) of Fig. 2, where 
the mass weighted stellar metallicity for the LAElg ranges 
between Z* = 0.04 — O.IZ©. Since the dust mass is calculated 
using the SFR averaged over the entire star formation his- 
tory of each progenitor {— Af,/i,), and since all the LAElg 
have ages between t, — 150 — 250 Myr as mentioned, progen- 
itors with the largest stellar mass are also the most dust rich 
with dust mass Mdust ~ 1O'''^~^'^M0, as seen from panel (c) 
of Fig. 2. To translate the dust mass into the color excess, 
we use the supernova extinction curve (Bianchi & Schneider 
2007). The color excess is then related to the escape fraction 
of continuum photons as E(B - V) = -2.51ogl0(/c)/11.08, 
using which E{B - V) ^ 0.03 - 0.048 for the LAElg 
as shown in panel (d) of Fig. 2; the results quoted here 
do not change using the Calzetti extinction curve where 
E{B -V) = -2.51ogl0(/c)/10.9. 



As for the IGM transmission for LAElg, it shows a 
scatter between Ta ~ 0.35 — 0.65 due to the varying con- 
tributions to the photoionization rates of the LAElg, from 
clustered sources. This clearly highlights the importance of 
such sources even for values of xhi as low as 10~^, as seen 
from panel (e) of Fig. 2; it is not necessarily the largest 
galaxies that have the largest contribution from clustering, 
as seen from the same panel. Finally, we show the observed 
Lya luminosity for the LAElg for xhi = 10~*: the scatter 
in Ta translates in a scatter into the observed Lya lumi- 
nosity, which ranges between 10*^~*^'^ergs~^ as shown in 
panel (f) of the same figure. As the value of xhi increases, 
both the value of Ta and La decrease. Therefore, the values 
shown in panels (e) and (f) represent the upper limits on 
the quantities shown. We can now compare the properties 
of the LAElg to those of the field LAEs in cosmological vol- 
umes {r^ lO^cMpc'^). The LAElg he at the lower mass end 
of the general LAE population, which have stellar masses 
~ igS-iO-^M© (panel al of Fig. 6, Dayal et al. 2009), stellar 
metallicity between 0.01 — O.5Z0 (panel a3 of Fig. 6, Dayal 
et al. 2009), dust mass between lO'^'^'^M© (Panel a of Fig. 
5, Dayal, Ferrara & Saro 2010) and an average color excess 
E{B - V) - 0.14 (Dayal, Ferrara & Saro 2010). 

Finally, we find that in our re-simulated volume (~ 
86cMpc^) the number density of LAElg is 10"°'^^ cMpc"^ 
for Xhi = 10"''. For a similar value of xhi, we sample 
a cosmological simulation snapshot at z '^ 6 (described 
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in Dayal, Ferrara & Saro 2010, which has a volume ol 
lO^cMpc'^), about ten thousand times by randomly placing 
the re-simulated volume within it. We find that the num- 
ber density of field LAEs is about two orders of magnitude 



lower, 



10' 



' cMpc . Even for xhi = 0.4, the number 



density of LAE^g (10"^''^ cMpc"^) is higher than that of 
the field LAEs by about one order of magnitude. Detection 
of such a large over-density of LAEs in a volume similar to 
the re-simulated one could be an excellent signature of a re- 
gion that could evolve to resemble the local group volume at 
z ~ 0. However, we caution the reader that to what extent 
such a result depends on the different resolutions of these 
two simulations still needs to be explored in full detail. 



5 CONCLUSIONS AND DISCUSSION 

We present the first work at high-z {z ~ 6) that couples state 
of the art high resolution gas-dynamical simulations of the 
local group run within the CLUES framework, with a LAE 
model to identify the local group progenitors (designated 
LAElg) that could be visible as LAEs at a time when the 
Universe was only about 1 Gyr old. The main results from 
this study are now summarized: 

• The first LAElg, which is a progenitor of the MW, 
appears at an average IGM neutral hydrogen fraction of 
Xhi = 0.4 (0.1) including (neglecting) the effects of clus- 
tered sources. As xhi decreases to 10~*, 5 progenitors each 
of the MW and M31 become visible as LAElg- 

• Ta shows a large scatter, ranging between 0.35 — 0.65 at 
Xhi ~ lO^'', due to the varying photoionization rate contri- 
butions from the clustered sources of the LAElg; clustering 
is imperative in making many progenitors visible in the Lya 
even for a fully ionized IGM. 

• The LAElg lie at the low mass end of the field LAEs: 
they have M. ~ IO^'^'^Mq, M^ust ~ 10''-^~^-^Mo and 
E{B -V) ^ 0.03 - 0.048. On the other hand, field LAEs 
are much larger and have M* ~ 1O'*"^°'^M0, M^ust ~ 
lO*-^-^M0 and an average E{B - V) = 0.14. 

• Finally, our results suggest an observable imprint of 
such high-z LAElg: the number density of LAElg is about 
two (one) orders of magnitude higher than that of field LAEs 
for Xhi = 10~^ (xhi ~ 0.4). Such high number- densities at 
2: ~ 6 would be excellent signatures of a region that could 
resemble the local group volume at 2: ~ 0. However, the 
extent to which such a result depends on the varying simu- 
lation resolutions still needs to be explored in more detail. 

The main caveats involved in this study are with regards 
to the distribution of dust (homogeneous/clumpy) in the in- 
terstellar medium of these high-z progenitors (e.g. Finkel- 
stein et al. 2009; Dayal, Ferrara & Saro 2010; Dayal, Maselli 

6 Ferrara 2011), and the effects of peculiar velocities (in- 
fiows/outflows into/from the galaxies) that can have enor- 
mous impact on Ta (Verhamme et al. 2006; Dayal, Maselli & 
Ferrara 2011). Uncertainties also remain on the value of the 
escape fraction of H i ionizing photons used. As discussed 
in the above mentioned works, it is unclear to what extent 
these effects can modify the observed Lya luminosity. It is 
hoped that a combination of the expected data from up- 
coming missions such as ALMA, MUSE and JWST will be 
fundamental in resolving such issues. 
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